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, In this article we derive formula for probability P(sup t<r (X(£) — ct) > 

u) where X = {X(t)} is a spectrally positive Levy process and c G M . 
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_cp_. 1 Introduction 

Levy processes arise in many areas of probability and play an important role 
among stochastic processes. Moreover they are a basic skeleton in many applica- 
tions in finance, insurance, queueing systems and engineering. Especially extreme 
values of stochastic processes are applicable in many fields. In this article we will 
consider spectrally positive Levy processes that is Levy processes with Levy mea- 
sure concentrated on (0, oo) . Let X = {X(t) : t > 0} be a Levy process with 
characteristic function of the form 



IE exp{iuX (t)} = exp < t 



lau + / (e iux - 1 - iuxl(x < 1)) Q(dx) 



(1) 

where a G R and Q is a Levy measure (this is a pure jump Levy process 
without Gaussian component). We assume that the Levy measure Q has a 
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density bounded on [e, oo) for every e > . We will investigate two cases that is 
the finite variation case and the compensated case. In the first case we assume 
Jo xQ(dx) < oo and the characteristic function of X we put 



TEi exp{iuX (t)} = exp < t 



{e iux - 1) Q{dx) 



(2) 



In the compensation case we have x Q(dx) = oo . 

Exact formulas for the supremum distribution of a stochastic process are 
known only in a few cases: Brownian motion, compound Poisson process, subor- 
dinators (see Takacs [II]), gamma process (see Dickson and Waters [7j), stable 
Levy process (see Bernyk et al. [3] and Michna [9]). In Takacs [H] a general- 
ization of the classical ballot theorem is provided to get the distribution of the 
supremum for the processes with interchangeable non-negative increments. 

In this article we will find the exact formula for the following probability 

T) = P(sup(X(t) - ct) > u) (3) 

t<T 

where T > , u > and c e M. . For c > probability ([3]) is the so-called 
finite time ruin probability. In the compensation case the drift c can be included 
in the parameter a . In the articles Dickson and Waters [7] and Michna [H] a 
special approximation of gamma process and a -stable Levy process by compound 
Poisson processes is used. In Michna [9] central limit theorem is applied to 
approximate a -stable Levy processes. The method to derive the formula for the 
supremum used in this paper is similar to the method applied in Dickson and 
Waters [7] and Michna [9] but here we apply a general approximation of a Levy 
process by compound Poisson processes which requires different arguments in the 
proof and get formulas for the probability (J3J) which in the compensation case 
seem to be a new result. 

We will use the results for the so-called finite time ruin probability for the 
classical risk process (the compound Poisson process). First let 

N(t) 

S(t) = Y,Y k -ct 

k=l 

be the classical claim surplus process, where {Yk\^ =l is an iid sequence of positive 
random variables, N(t) is a Poisson process and c > . The first formula below 
is going back to Cramer (see e.g. Asmussen and AlbrecherpQ) and the second one 
is the so-called Seal's formula (see Seal [13], originating from Prabhu [ID])- Let 
ip(u, T) = P(sup t<T S(t) > u) . 

Theorem 1. For T > , u > and c > we have 
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(V 



1-<0(O,T) = P(supS(*) = 0) = — / P ( > ^ Y k < x | dx 



1 r cT ( N{T) \ 



(a) 



1 - ^(m,T) = P(sup£(i) < u 

t<T 



= P(S(T) < u) - c [ (1 — ip(0, T — s))lP(S(s) G [u,u + ds]) 
Jo 

= X>(S(T) <u)-c I (1 - V>(0, T - s))/(w + cs, s) rfs , 
Jo 

where f(x, s) is the density function of Ylk=i ^ > provided that it exists. 

Remark 1. A compound Poisson distribution has an atom at thus the above 
function f is a density of the positive part of the distribution support. 

We will apply approximation of a Levy process X by compound Poisson 
processes. Let us put 

N e (t) = ^2aX(s)1(AX(s) > e) . 



s<t 



The process N e is a compound Poisson process with positive jumps. 

Proposition 1. Let X be a spectrally positive Levy process with a Levy measure 
Q . Then 

(i) for xQ(dx) < oo 

N e (t)^X(t); (4) 

(ii) for Jq 1 xQ(dx) = oo 

N e (t) + (a- £ xQ(dx))t^X(t), (5) 
as e I a.s. in the uniform topology. 

Proof. The assertions follow from Levy-Ito representation see e.g. Sato [12J. □ 

2 Main result 

Now we use the above results to derive a formula for probability (j3]) (for a similar 
treatment for gamma process see Dickson and Waters [7] and for a -stable case 
Michna [9]). 
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Theorem 2. Let X be a spectrally positive Levy process, T > and u > then 
(i) for xQ(dx) < oo and c > 

i r T 

(0, r) = P(sup(X(t) - ct) > 0) = 1 - — / P X(T) < x) 



and 



*(u,T) = P(sup(X(t) - ct) > u) 

t<T 

= P(X(T) -cT > u) + 

T f(u + cs,s^ /' C(T - 1 



T-s 



ds / P(X(T- s) < x)ofo, 



'0 J ° jo 

where f(x,s) is the density function of X(s) , provided that it exists; 

(ii) for xQ(dx) = oo and a6l 



*(0,T) = P(supX(t) > 0) = 1 



and 



*(u,T) = P(supX(t)>u) 
= P(X(T) > it) + 

/ ^^-ds [ 1P(X(T-s)<x)dx, (6) 

JO -L — S J -oo 

where f(x, s) is the density function of X(s) , provided that it exists and 
the cumulative distribution of sup t<T X(t) is continuous. 

Remark 2. The formula in (i) is similar to the formula of Takdcs [Lffl valid for 
the processes with interchangeable non-negative increments. 

Remark 3. The formula for ty(u,T) can be unified for the both cases because in 
the first case the support of the distribution of X(t) is in [0, oo) and the formula 
(0|) is valid for the case (i) too. 

Remark 4. Formula (Q|) is valid for Wiener process too (see Michna JM/). 

Remark 5. In the case (i) the distribution of sup t<T (X(t) — ct) has an atom at 
u = . 

Remark 6. The case (i) for c < is trivial. 
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Proof. First we show the assertions of (i). Let {e ra }^ =1 be a decreasing sequence 
of positive numbers such that lim^oo e n = . Then it is easy to notice that 
N en (t) < AT (t) for all t > and n. Thus 



and 



sup(iV en (t) - ct) < sup(iV en+1 (t) - ct) 

t<T t<T 



lim sup(iV £n (t) - ct) = sup(X(i) - ct) . 



Hence for all u6l we have 

lim P(sup(iV en (t) - ct) <u) = P(sup(X(t) - ct) < u) . 

Using (i) of Theorem [TJ we get for c > that 

l-tp e (0,T) = P(sup(iV £ (t) -ct) < 0) 

t<T 

= ^ P (N e (T) < x) dx (7) 



1 /- cT 
cT Jo 



where in the last line e 4- . Thus we get the first assertion of (i). 

Now let us consider (ii) of Theorem[T]for the process N e . The random variable 
N e (t) has a density by the assumption that the Levy measure Q has a bounded 
density on [e, oo) . Thus we have for c > and u > 

1 - ip e (u,T) = P(sup(iV £ (f) - ct) < u) 

t<T 

= P(iV e (T) - cT <u) - c I {l-i/j e (0,T-s))'P(N e {s)-cse[u,u + ds)) 

Jo 

= P(iV e (T) -cT <u)- 

nT 1 MT-s) 

/ P(iV e (s) -cs e [u,u + ds]) / P (iV e (T — s) < x) dx 

Jo T — s J 

-»■ P(X(T)-cT<u)- 

/ P(X(s) -cs G [u,u + ds]) / JP(X{T-s) <x) dx, 

Jo T — s J 

where in the second last equality we substituted (JZJ) and in the last one e | . 
This way we obtain the second assertion of (i). 

Now let us prove the assertions of (ii). We put c e = —a + xQ(dx) and 
notice that lim^o c e = oo . By the property (ii) of Proposition CE] we have 

limsup(./V € (i) — c t t) = supX(t) 

40 t<T t<T 
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which implies that 

limP(sup(iV e (t) - c e t) < u) = P(supX(t) < u) 

40 t<T t<T 

for every u G M where P(sup i<r X(t) < u) is a continuous function of u . Thus 
for a sufficiently small e we have c e > and we get 

1 - V^(0, T) = P(sup(iV £ (t) - c e t) < 0) 

t<T 



-±= [ E P (N e (T) < x) dx 

CeT Jo 
1 f° 

— / P (N e (T) - c e T < x) dx (8) 



-)- 



where in the last step e 4- and the last integral tends to 

J"° P (X(T) < x) dx which is finite by the fact the Levy measure Q is concen- 
trated on (0, 00) . For u > we obtain 

1 - tf) € (u, T) = P(sup(iV £ (t) - c e t) < u) 

t<T 

= P(JV e (T) - c e T < u) - c e [ {l-iJj e (0,T-s))P{N e (s)-c e se[u,u + ds)) 

Jo 

= P(iV e (T) - c t T <u)- 

[ 7 ^—P(N t (s) -c e s e [u,u + ds]) ■ 
Jo T-s 

[ P (N e (T — s) — c e (T — s) <x) dx 

J-cAT-s) 



-Ce(T-s) 

->■ P(X(T) < U 



f T ds r P (X(T -s)<x) dx 

Jo T-s J^ 



where in the second last equality we substituted © and in the last one e 4- . □ 

We can give the following sufficient condition for the Levy process X to have 
the cumulative distribution of sup t<T (X (t) — ct) continuous. 

Proposition 2. If X is a Levy process such that 

X(t) 

limsup^^>0 (9) 

40 t K 

a.s. for k < 1 then the cumulative distribution of sup t<T (X (t) — ct) is continu- 
ous. 

Remark 7. Conditions for the property (Tj|) are given in Bertoin et al. J^j. 
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Proof. By ([9]) we get 

X(t)-ct 
hm sup > (J 

a.s. which gives that the process X(t) — ct is above the barrier u = infinitely 
many times in every neighborhood of t = a.s. Thus 

P(sup(X(t) - ct) = 0) = . (10) 

For «>0 we define the stopping time t(w) = inf{£ > : X(t) > u] and by the 
strong Markov property we get that the process Y(t) = X(t + t) — X{j) is a 
Levy process with the same distribution as the process X on the set {r < oo} . 
Thus using ( flOj) we obtain that 

P(sup(X(t) - ct) = u) = . 

t<T 

□ 

For stable Levy processes the continuity of the cumulative distribution of 
sup t<T (X(£) — ct) can be shown using the law of the iterated logarithm for small 
times (see Zolotarev [15]). Moreover Theorem |2] can be applied to many other 
Levy processes such as the spectrally positive generalized hyperbolic processes 
(e.g. the inverse Gaussian Levy process) introduced by Barndorff-Nielsen and 
Halgareen [2], the spectrally positive Levy processes proposed by Carr et al. 
[6] or tempered stable processes introduced by Koponen [8] and generalized by 
Rosihski [TT]. Recently the inverse Gaussian Levy process is extensively used 
in financial modeling. Thus let us derive the exact formula for the supremum 
distribution of the inverse Gaussian process. 

Example 1. Let us consider the inverse Gaussian process X(t) with the following 
density 



[ 2 

where 5 > and 7 > 0. The process X is a subordinator. By Theorem^ (i) 
we have for c > 



P(sup(X(t) - ct) > 0) 

t<T 

-cT nX 



l ~ 5 ^m[ dx J* y~ 3/2 {-l^y" 1 + ^} d y 

J^ T (cTx^ 2 - x- 1 / 2 ) exp j-i^TV 1 + 7^)} dx , 



5el 8T rcT 

1 - 



cV2n 
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where in the last step we integrate by parts. For u > we consider the second 
assertion of Theorem^ (i) and similarly we obtain 



P(sup(X(t) - ct) > u 

t<T 



+ 



6 2 e yST f T 



r x- 3 / 2 exp(-J(5 2 T 2 x- 1 + 7 2 x)) dx 

Ju+cT 12 J 

J s(u + cs)' 3/2 exp j-^(<5V(u + cs)' 1 + 7 2 (m + cs))J ds 



2n 

• J (c(T - s)x-' i/2 - x~ 1/2 ) exp |-- (5 2 (T - s) 2 ^- 1 + 7 2 x)| dx . 

One can investigate a tempered stable Levy process X with the Levy measure 

Ae~ Xx 

Q(dx) = — I{x > 0} dx , 

where A > , A > and < a < 2 . The closed form of the density function 
of X(t) is not known but the characteristic function of the distribution of X(t) 
can be given in a closed form see Carr et al. [5]. Thus we can calculate the 
inverse Fourier transform numerically to get the density function of X(t) . For 
< a < 1 we apply Theorem [2] (i) and for 1 < a < 2 we use Theorem |2] (ii) to 
get the formula for the supremum distribution of the process X . 
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